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Summary 

The basic functions embedded within a user- 
friendly software package (MATRIXx) are used to 
provide a high-level systems approach to the analy- 
sis of linear control systems. Various control-system 
analysis configurations are assembled automatically 
to minimize the amount of work by the user. 
Interactive decision making is incorporated via menu 
options and at selected points, such as in the plot- 
ting sections, by inputting data. One section in this 
paper describes the frequency and steady-state anal- 
ysis. There are five evaluations — the singular-value 
robustness test, the singular-value loop-transfer fre- 
quency response, the Bode frequency response, the 
steady-state covariance analysis, and the closed-loop 
eigenvalue calculations. Another section describes 
time-response simulations. These simulations, which 
are linear except for actuator saturation, may have 
either a continuous plant or a discrete plant and in- 
clude a generic controller structure for both feedback 
and feedforward dynamics. A time response for a 
random white-noise disturbance is available. This 
paper describes the configurations and key equations 
used for each type of analysis, the restrictions that 
apply, the type of data required, and a sample prob- 
lem. One approach to integrating the design and 
analysis tools is also presented. 

1.0. Introduction 

Prior to the early 1960’s, hand calculations using 
slide rules and spirules were the major tools for con- 
trol design; analog computers were the major tools 
for simulation. During the 1960’s and 1970’s, as 
control design and analysis techniques became more 
complex, use of more powerful digital computers be- 
came necessary. Specialized control-system design 
and analysis programs, usually written in Fortran, 
were developed for use on digital computers. These 
programs were significantly more powerful than pre- 
viously available techniques, but quite often the 
designer spent a large amount of time debugging soft- 
ware rather than solving the control problem. Typ- 
ically, the control designer needed a programmer to 
interface with the computer because of the program- 
ming complexity. 

Computer-aided control-system design and analy- 
sis tools have emerged during the past decade. These 
programs provide a new generation of user-friendly 
tools that allow the designer to spend more time 
on the primary task of control-law design. A sum- 
mary of 20 computer-aided control-system design 
software packages is shown in reference 1, and ad- 
ditional packages have been developed since this ref- 
erence was published. A user-friendly software pack- 


age described in the software summary MATRIXx 
(ref. 2) has been purchased by various organizations 
at Langley Research Center. This paper is not an 
endorsement of MATRIXx over any other software 
package. However, since MATRIXx is the software 
package used by our organization, the control-system 
analysis tools described in this paper use the funda- 
mental features of this computer program. 

The program MATRIXx has the basic func- 
tions that allow implementation of a high-level sys- 
tems approach to control-system analysis. These 
basic functions include matrix algebra, transfor- 
mations between frequency and time domain, fre- 
quency responses, time simulation, plotting fea- 
tures, and structured programming capability. Also, 
SYSTEM-BUILD (an optional software module, for 
use with MATRIXx, that is referred to as BUILD 
in this paper) provides an environment for build- 
ing, modifying, and editing simulation models. Ver- 
sion 6.0 of BUILD has a programming capability that 
allows building blocks to be automatically assembled 
into various system configurations. 

The approach described in this paper is to use 
the basic functions within MATRIXx to provide a 
systems approach to the analysis of linear control 
systems. Using this approach, designers can perform 
various analyses quickly and accurately. The idea is 
to minimize the calculations and matrix manipula- 
tions by the designer and still provide interactive flex- 
ibility. Decisions are made via menu-driven options, 
which have been programmed into the major parts 
of the software. At selected points in the software, 
mainly in the plotting sections, the designer can also 
interact with the control-system analysis tools by in- 
putting data. For example, the designer may input 
data to create a vector of variable-density-frequency 
points for selected frequency-response evaluations. 

All the evaluations are developed for discrete con- 
trol systems. The main input data required are state- 
space quadruples in discrete form for both a plant 
and a controller. However, continuous systems can 
also be evaluated, as is described in this paper. The 
program BUILD is used to construct the configura- 
tions needed by the frequency and time analyses. 
The system quadruple generated within BUILD is 
then used in the appropriate evaluation. 

Sections 3.0 and 4.0 of this paper describe the 
control-system analysis capabilities that have been 
developed. Frequency and steady-state analyses are 
described in section 3.0. Five evaluations presently 
available are the singular-value robustness test, the 
singular- value loop-transfer frequency response, the 
Bode frequency response, the steady-state covariance 
analysis to an external disturbance, and the closed- 
loop eigenvalue calculations. The discussion of the 


five evaluations includes a description of menu op- 
tions, key equations, and control-system analysis di- 
agrams. There are also subsections that describe 
the plotting of the frequency data and describe an 
optional user file. 

The description of time analyses is in section 4.0 
of this paper. Time-response simulations are linear 
except for actuator saturation, which can be included 
at the option of the user. The user can also choose 
between a continuous plant and a discrete plant 
in the simulation. Subsections describe a generic 
controller structure that allows feedback and feed- 
forward dynamic controllers, signal-generator op- 
tions, and control-loop construction. Other sub- 
sections describe the time response to external 
disturbances and the plotting capabilities for the 
time-response data. 

An approach that allows integration of a control- 
system design program and the analysis tools is 
described in section 5.0 of this paper. This approach 
can be used to reduce the total overall design time 
by allowing the designer quick and easy access to the 
analysis tools. 

A sample problem is described in section 6.0 
of this paper. Details of the sample problem are 
shown in the appendixes; appendix A contains the 
input data and displayed output of the frequency- 
analysis evaluations, and appendix B contains similar 
contents for the time-analysis evaluations. 

One objective of this development is to enable 
the user to assemble and analyze various control- 
system configurations using existing MATRIXx func- 
tions, user-defined functions (UDF), command files, 
and macros. A goal is to have, in a single software 
package, user-friendly control-system analysis tools 
that perform a variety of analyses. These analysis 
tools are easily integrated with control-law synthesis 
tools to form a complete environment. Within this 
paper, examples of flight control problems are used, 
but the approach is generic and should be useful to 
other control disciplines. All the control-system con- 
figurations assume positive feedback with the nega- 
tive sign included within the control block. Features 
described are applicable to MATRIXx version 6.0, 
although upgraded versions will probably necessitate 
some changes in the software; the basic approach 
should remain applicable. All UDF’s end with the 
notation “.FUN”, all command files end with the no- 
tation “.CMD”, and all plot command files end with 
the notation “.PLT”. 

2.0. Nomenclature 

A state matrix for continuous 

system 


(A,B,C,D) 

quadruple matrices for contin- 
uous system 

ANSW 

variable to differentiate 
whether same or different sig- 
nals are used for all models 


transfer-function matrix 

B b ,B f 

control matrices (for feedback 
and feedforward signals, 
respectively) used in generic 
controller 

C 

state-space quadruple matrix 
of controller or state output 
matrix of a continuous system 

C1,C2, C3 

matrices used in generic 
controller to describe different 
blocks 

D(,,Dy 

input-to-output transfer ma- 
trices (for feedback and feed- 
forward signals, respectively) 
used in generic controller 

D, 

diagonal scaling matrix 

DBODE 

MATRIXx discrete Bode 
response function 

DELT 

discrete sampling period, sec 

DELTIME 

time step for each element in 
output time vector, sec 

DFP 

D avidon- F le t cher- Powell 

DIMS 

row vector containing scalar 
input data 

DLSIM 

MATRIXx function for dis- 
crete time response to general 
inputs 

E 

input-to-output transfer 
matrix for discrete system 

EVALCL 

table containing vectors of 
closed-loop eigenvalues and 
damping ratios 

(F,G,H,E) 

quadruple matrices for a 
discrete system 

F 

state transition matrix of a 
discrete controller 

FREQ 

MATRIXx frequency-response 
function 

G 

control matrix for a discrete 
system 
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GAINM 

gain margin, dB 

GUSTAT 

statistical data (standard 
deviation and mean) generated 
from time-simulation data 
with random input on first 
external disturbance path 

GWSTAT 

statistical data (standard 
deviation and mean) generated 
from time-simulation data 
with random input on second 
external disturbance path 

H 

state output matrix for a 
discrete system 

I 

identity matrix 

inf 

infimum (greatest lower 
bound) 

j 

= V=i 

LBOUND, 

UBOUND 

lower and upper bounds for 
actuator saturation 

NMOD 

number of models 

NS 

number of states for SDP or 
SCP 

NSCP 

number of states for SCP 

NSCl, NSC2, 
NSC3 

number of states for SCI, SC2, 
and SC3, respectively 

NSDC 

number of states for SDC 

NSGDP 

number of states for SGDP 

NU 

number of control inputs (used 
in Bode response) 

P 

state-space quadruple matrices 
of plant 

PHASEM 

phase margin, deg 

Q 

weighting matrix 

RMS 

root mean square 

ROUND 

MATRIXx function for com- 
puting an integer closest to 
input number 

SCI, SC3 

feedforward quadruples for 
generic controller 

SC2 

feedback quadruple for generic 
controller 

SCP 

quadruple for continuous plant 
model 


SDC 

quadruple for discrete con- 
troller model 

SDP 

quadruple for discrete plant 
model 

SGDP 

quadruple for a discrete 
model made up of a combined 
airplane and disturbance 
model 

SIGI 

minimum singular value of 
inverse return difference 
matrix at plant input 

SIGO 

minimum singular value of 
inverse return difference 
matrix at plant output 

SISO 

single-input, single-output 

SVLT 

singular- value loop transfer 

SVLTI 

singular value of loop transfer 
at plant input 

SVLTO 

singular value of loop transfer 
at plant output 

T 

discrete sampling period, sec 

TMAX 

maximum time of simulation 
run, sec 

T S 

output time step 

U 

input disturbance 

UDF 

user-defined function 

X 

steady-state covariance matrix 

XURMS 

RMS values for states with 
random disturbance on first 
external disturbance path 

XWRMS 

RMS values for states with 
random disturbance on second 
external disturbance path 

Y 

output 

YURMS 

RMS values for outputs with 
random disturbance on first 
external disturbance path 

YVV 

row vector that defines outputs 
for feedback to controller 

YWRMS 

RMS values for outputs with 
random disturbance on second 
external disturbance path 

z 

discrete frequency parameter 


A stable transfer matrix that 

is either fully populated or 
diagonal 

S srn stability margin 

7] random white-noise input 

ji structured singular value 

a maximum singular value 

uj angular frequency, rad/sec 

u;g US t colored noise output from gust 

spectrum (fig. 8) 

Superscripts: 

T transpose 

— 1 inverse 

Special notation: 

.CMD suffix for command files 

.FUN suffix for user-defined 

functions 

.PLT suffix for plot files 

1 1 S 1 1 2 2- norm of a matrix S 

3.0. Frequency Analyses 

The types of analyses described in this section are 
frequency responses, for several control-system con- 
figurations, and steady-state solutions. Each control- 
system configuration is composed of linear models 
that represent a discrete controller and a discrete 
plant; these models are assembled by using the pro- 
gramming capability within BUILD. All the analy- 
sis routines have been developed to handle several 
systems automatically in one computer run, which 
allows comparison of different control systems at dif- 
ferent flight conditions. One restriction is that all 
systems being compared must have the same num- 
ber of states, inputs, and outputs. However, models 
with different numbers of states can always be com- 
pared by adding dummy states with stable eigen- 
values to the lower order models so that the order 
of each model is the same. 

Continuous-system analyses can also be accom- 
plished by using these control-system analysis tools. 
A continuous system can always be discretized at a 
high frequency compared with the bandwidth of in- 
terest and then analyzed as a discrete system. A 
Tustin transformation (ref. 3) is a good discretiza- 
tion approach, because it does not add the phase lag 
associated with a zero-order hold. TUSTIN. FUN has 
been developed to allow users to analyze continuous 


control systems. It is suggested that the sampling fre- 
quency be at least four times greater than the highest 
frequency of interest. 

The structure for the frequency- analysis sec- 
tion is shown in figure 1. FREQANAL.CMD 
serves as an executive and calls other command 
files (FREQEVAL.CMD, PLOT.CMD, and 
USERFILE.CMD) that execute each of the main 
parts of the frequency- analysis software. All the con- 
figurations for the frequency-evaluation analyses are 
assembled and evaluated using FREQEVAL.CMD. 
Data that are generated can then be displayed using 
PLOT.CMD. USERFILE.CMD is an optional user- 
developed file that can be used to create all the data 
needed for the frequency-analysis and plotting com- 
mand files. When USERFILE.CMD is selected via 
a menu option, the user is prompted to specify the 
name of the desired file. As an example, the first au- 
thor has developed a command file (SYSTEM.CMD) 
that uses the input and output data from a control 
design program and generates discrete quadruples of 
the plant and controller. Whether a user file is used 
or not, any of the main command files can be run 
separately if the required data are available. 

The following subsections describe the contents of 
each command file and associated UDF. 

3.1. Frequency Evaluation 

Five different frequency and steady-state analyses 
are presently available: 

1. Singular- value robustness test 

2. Singular- value loop transfer 

3. Bode frequency response 

4. Steady-state covariance analysis to external 
disturbance 

5. Closed-loop eigenvalues 

The command file FREQEVAL.CMD manages the 
various analyses. Figure 2, the frequency evaluation 
flowchart, shows the UDF’s that have been developed 
and the flow structure. There are many menus in 
the software. A typical frequency-evaluation menu 
(fig. 3) is invoked in FREQEVAL.CMD to select any 
of the five analyses. 

A menu option is also available to change the fre- 
quency range, since it is typical to perform some 
of the frequency- response tests over different fre- 
quency ranges. NOMEG. FUN, which is used for 
all the frequency-response tests, has been developed 
to create a vector of frequency data points that are 
linearly spaced on a logarithmic scale. The inputs to 
NOMEG. FUN allow multiple frequency ranges with 
different data-point densities for each range, speci- 
fied as the number of points per decade. The vector 
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of frequency data points is then composed of sub- 
vectors, each with its specified density of frequency 
points. This feature is important when evaluating 
a control system with a lightly damped eigenvalue, 
because additional data points can be calculated in 
the frequency range where the control-response vari- 
able is varying relatively fast. In contrast, fewer data 
points are needed in frequency ranges with slowly 
varying control responses. This flexibility results in 
reduced time to perform the analysis and still allows 
for the important response characteristics. 

The program BUILD is used to assemble the sys- 
tem configurations in all the evaluations described. 
The programming capability of BUILD allows an au- 
tomated approach to the assembly of the system con- 
figuration, which speeds the evaluation process and 
helps minimize errors. The only data required are 
the discrete system quadruple matrices (F,G,H,E) 
for both the controller and plant and the number 
of states in each quadruple. The only restrictions 
are that the number of plant outputs must equal 
the number of controller inputs and the number of 
controller outputs must equal the number of plant 
inputs. The program can calculate the number of in- 
puts and outputs to each dynamic system and make 
the appropriate connections. There is no limit to the 
number of systems that can be evaluated in one com- 
puter run. The five types of analysis are described in 
more detail in the subsections that follow. 


The state-space quadruple matrices of the controller 
and plant are represented by C and P, respectively, 
and (i represents the structured singular value. 

The expressions shown in equations (2) and (3) 
are for a closed-loop system with input disturbance 
U and output Y. These expressions are implemented 
using the input-to-output map shown in figure 5. 
SV.FUN is supplied with quadruples of the discrete 
plant and discrete controller along with the sam- 
pling period and outputs closed-loop-system quadru- 
ple matrices that represent the transfer function for 
the configuration shown in figure 5. SV.FUN uses 
BUILD to create this configuration. After the first 
model is complete, the procedure is repeated for 
other models, except that BUILD is used only to up- 
date the quadruples. SSVIO.FUN is supplied with 
the quadruples generated in SV.FUN and outputs 
the singular values of the system being processed 
along with the minimum singular value of the inverse- 
return difference matrix. Within SSVIO.FUN, the 
closed-loop quadruple (F,G,H,E) representing the 
transfer function is then split into components to 
compute the complex frequency response for A s in 
the discrete domain (ref. 3) as follows: 


A s = H(zl - F) -1 'G + E (4) 


3.1.1. Singular-value robustness test. The 
option exists to calculate the stability margins of a 
closed-loop system for a multiplicative error at either 
the plant input or plant output. The results of this 
test give an indication of system stability when all 
loops are opened with simultaneous gain and phase 
changes in all loops. An example of the multiplicative 
error at the plant input is shown in figure 4, where 
a stable transfer matrix A is assumed to be acting 
simultaneously on all input paths. The closed-loop 
system is guaranteed to be stable if (ref. 4) 

d“(A) < 6 S m (I) 


where a(A) is the maximum singular value and 6 sm is 
the stability margin. The calculation for 6 sm depends 
on whether A is fully populated or diagonal (ref. 5) 
as follows: 


l 

_ ^[(l-CP)- 1 CP] 

1 

[(i — cp) _ 1 cp] 


(for A fully populated) (2) 

(for A diagonal) (3) 


z = exp(juT) (5) 


where w is the angular frequency, T is the discrete 
sampling period, and I is an identity matrix. A 
menu option is available to allow the solution of equa- 
tion (4) directly for the unstructured singular value 
shown in equation (2) or to allow the computation of 
the structured singular value shown in equation (3). 
One shortcoming of the unstructured-singular-value 
test is the conservativeness of the results when con- 
sidering multiloop-system stability margins (ref. 6). 
The singular value changes with system units (de- 
grees, radians, etc.). This scaling problem contrasts 
with the classical Bode test for stability, which is in- 
variant with changes in the system units. Reference 5 
shows that a scaling matrix can be used to eliminate 
the scaling problem. The solution for /x requires the 
calculation of a diagonal scaling matrix D s at each 
frequency as follows: 


m(F,G,H,E) < inf^D.A.DJ 1 ) (6) 

Ds 
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Unfortunately, D 5 must be solved iteratively 1 as an 
optimization problem by minimizing a with respect 
to the components of T> s . When the number of paths 
is three or less, relation (6) becomes an equality, but 
when the number of paths is greater than three, the 
upper bound for a becomes a reasonable approxima- 
tion to fjL. Two convergence algorithms are available; 
the default algorithm is a Davidon-Fletcher-Powell 
(DFP) algorithm (ref. 7), and the other is a modified 
Newton-Raphson algorithm. Based on the authors’ 
experience, the DFP operates significantly faster and 
has been successfully used on several different prob- 
lems. The second algorithm requires the calcula- 
tion of a Hessian matrix, which takes a considerable 
amount of time. SVDFP.FUN and SVGNR.FUN 
have been developed to execute the DFP and modi- 
fied Newton-Raphson algorithms, respectively. 

A menu option allows selection of the iteration 
rate for calculating D 5 . The most accurate choice is 
to compute D s at every frequency. To reduce cal- 
culation time, the user can compute D s less often. 
For example, if the default value of five is chosen 
for the D s computation, the same D 5 is used for 
five consecutive frequencies. The five frequencies in- 
clude the two frequency points immediately preced- 
ing and immediately after the frequency at which the 
computation takes place. 

3.1.2. Singular-value loop transfer. The 

singular-value loop transfer (SVLT) is an open-loop 
frequency-response test for which all loops are bro- 
ken either at the plant input or the plant output. 
This test gives an indication of open-loop response, 
such as the crossover frequency and high-frequency 
attenuation. In a single-input, single-output (SISO) 
system, the SVLT is identical to the magnitude re- 
sponse in the Bode test. SVLT. FUN has been devel- 
oped for the SVLT test. Within SVLT.FUN, BUILD 
is used to assemble two configurations (fig. 6), one 
for the loops broken at the plant input and the other 
for the loops broken at the plant output. Both the 
input and output SVLT responses are automatically 
generated over the frequency range of interest, since 
the calculation is relatively fast. The configuration 
for the plant input is assembled and linearized, and a 
singular-value frequency response is then computed 
using the MATRIXx function FREQ. The procedure 
is repeated for the output configuration. After the 
first model has been analyzed, the procedure is re- 
peated for other models; BUILD is used only to up- 
date the plant and controller matrices. 


1 Work done by Barton J. Bacon of Langley Research Center 
as part of a Ph.D. thesis entitled “Order Reduction for Closed- 
Loop System” for the School of Aeronautical and Astronautical 
Engineering, Purdue University. 


3.1.3. Bode frequency response. The Bode 
frequency response is a classical stability test for 
SISO dynamic systems. BFR.FUN is supplied with 
system quadruples for a discrete plant and discrete 
controller, with a frequency vector, and with a dis- 
crete sampling time. The outputs of BFR.FUN are 
vectors of gain and phase and data for gain and phase 
margins. In BFR.FUN, BUILD is used to construct 
the various configurations; each plant input is opened 
consecutively with all other plant inputs closed. This 
construction is illustrated in figure 7, where one loop 
is shown broken at the plant input and the other 
loops remain closed. The box with 1” in figure 7 
is included to create a 180° phase shift that normally 
occurs with negative-feedback control systems. In 
this paper, positive feedback is assumed. 

A system quadruple is formed and a frequency 
response is computed using the MATRIXx function 
DBODE for each opened single-loop configuration. 
In the example of figure 7, NU frequency responses 
are computed for each model evaluated. After one 
model has been analyzed, the procedure is repeated 
for other models. 

DBODE in MATRIXx version 6.0 allows the 
phase angle to exceed the principal angle range of 
360° and results in decreased resolution on phase- 
angle plots. PRA.FUN has been developed to 
convert the phase-angle vector data to a princi- 
pal angle between 0° and —360°. A second UDF, 
MMARGIN.FUN, accepts the principal phase angle 
from PRA.FUN, along with amplitude data from 
DBODE, and then calculates the gain and phase 
margins. Linear interpolation is used to calculate the 
amplitude crossover frequency and the correspond- 
ing phase angle, to which 180° is added to calculate 
phase margin. Linear interpolation is also used to 
calculate the frequency at which the phase angle be- 
comes — 180°. The gain margin is then calculated 
from the amplitude data at the corresponding fre- 
quency. A stable control system must have a positive 
phase margin and a negative gain margin (in dB). 

3.1.4. Steady-state covariance analysis. 

Typically, a control engineer is interested in eval- 
uating the response of a closed-loop control sys- 
tem to random external disturbances such as wind 
gusts. SSCOV.FUN calculates the state and out- 
put steady-state covariances for two random exter- 
nal disturbances. The data required for SSCOV.FUN 
are the sampling period and discrete-system quadru- 
ple matrices. These matrices represent the con- 
troller and a combined plant composed of the distur- 
bance model and original plant. The outputs from 
SSCOV.FUN are the state and output RMS values 
that result from the two external-disturbance inputs. 
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A representative block diagram with a gust spec- 
trum as the external-disturbance model is shown in 
figure 8. 

Input signals to the combined plant are the con- 
troller feedback commands plus two signals for the 
random external disturbances; the output signals 
from the combined plant are the measurements going 
to the controller plus two outputs from the distur- 
bance model. Input and output restrictions are that 
the two random external-disturbance signals must be 
the last two input paths going to the combined plant 
and that the outputs from the disturbance model 
must be the last two outputs from the combined 
plant. The reason for this restriction is that all con- 
nections, except the last two inputs and outputs, are 
made between the plant and controller. Furthermore, 
the disturbance-model outputs are later normalized 
to unity. If there is only one input disturbance of 
interest, the second disturbance path can always be 
used as a dummy input. 

BUILD is used to construct the configuration of 
figure 8, which is then linearized to form the discrete- 
system quadruple (F,G,H,E). A MATRIXx function 
is used as follows to solve the discrete Lyapunov 
equation: 

X = FXF r + GQG r (7) 

where F is the discrete state transition matrix, G 
is the discrete control matrix for the random inputs, 
X is the steady-state covariance matrix, and Q is a 
weighting matrix. The state vectors that are output 
from SSCOV.FUN are composed of the square root 
of the diagonal elements of X. The output vectors are 
calculated by taking the square root of the diagonal 
elements of HXH , where H is the output matrix 
and E is assumed to be zero. The steady-state 
covariance is calculated separately for each input, 
and the RMS data are then normalized by the RMS 
values of the external-disturbance-model outputs. 

3.1.5. Closed-loop eigenvalues. A capabil- 
ity to calculate the closed-loop eigenvalues from the 
quadruple matrices of the discrete plant and discrete 
controller has been developed using CL. FUN. This 
function is used mainly as a check for the user, who 
typically calculates the closed-loop eigenvalues dur- 
ing the design process. Using BUILD, the closed- 
loop configuration is constructed and then linearized 
to form a quadruple matrix. The closed-loop discrete 
eigenvalues are calculated along with the equivalent 
continuous eigenvalues and damping ratios. This lat- 
ter step is accomplished in ZTOS.FUN. Vectors rep- 
resenting the discrete eigenvalues, the corresponding 
continuous eigenvalues, and the damping ratios are 
displayed in a table. 


3.2. Plotting 

Four of the five analyses described in section 3.1 
(Frequency Evaluation), with the closed-loop eigen- 
value calculation being the exception, have plotting 
command files that can be executed. PLOT.CMD 
(fig. 1) serves as the executive for the specialized plot 
command files. Several menus are defined and in- 
voked within PLOT.CMD. One menu is for selecting 
the specific plot command file to be executed. The 
choices are as follows: 

1. SV.PLT 

2. SVLT.PLT 

3. BFR.PLT 

4. SSCOV.PLT 

These plotting files are used to view the results 
from the corresponding analysis (items 2 through 
5, respectively, in fig. 3) selected in the frequency- 
evaluation menu. A second menu allows the selec- 
tion of six different screen options. A third menu 
allows hard-copy plots to be made. The specific plot 
command files are described in detail below. 

The minimum singular value of the inverse return 
difference matrix can be displayed on a log-log plot 
using SV.PLT. A menu option is invoked in SV.PLT 
to allow the user to display the singular values for 
either the plant input or the plant output. The 
available selections are determined by the type of 
analysis that has been performed. Multiple-model 
data can be put on the same plots when comparing 
various models. 

SVLT.PLT allows the user to plot the singular 
values of the loop-transfer function. A menu option 
allows the user to display the singular values for 
either the plant input or the plant output. The 
user has the option of displaying multiple-model data 
on a common graph. The user also may choose 
the singular values that are to be displayed on this 
common graph. All the singular values are displayed 
on a log- log plot. 

Bode frequency-response data are displayed on 
semilog plots by using BFR.PLT. A menu option 
allows the user to display either the amplitude in 
decibels, the phase angle (deg), or a Nyquist polar 
plot. The magnitude on the Nyquist plot is in actual 
units, as opposed to the desired units of dB. A dB 
scale is not used because MATRIXx polar plots do 
not allow negative numbers. Multiple models can be 
displayed on the same plot, and the desired open-loop 
data can be selected and displayed. 

The steady-state covariance data generated in 
SSCOV.FUN are displayed by SSCOV.PLT by using 
bar charts. A menu option allows the user to select 
plots of either the states or output variables for one 
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of the two external-disturbance paths. Two different 
prompts allow the user to enter, in vector form, the 
models to be used for the plot and to select the vector 
of elements (states or outputs) that will be displayed. 

3.3. User File 

One menu option in FREQANAL.CMD (fig. 1) 
is the choice of a user file, which is an optional 
file to be developed by the user. When this option 
is selected, the user is prompted to enter a com- 
mand file name. The user- file option has been in- 
cluded mainly as a convenience to allow the user a 
method of automatically preprocessing data for use 
in the frequency-evaluation file. This option does 
not have to be invoked if the user has already loaded 
the required data into MATRIXx prior to executing 
FREQANAL.CMD. (See appendix A and “Sample 
Problem,” section 6.0.) 

SYSTEM.CMD is an example of how the user- 
file option can be utilized. SYSTEM.CMD was de- 
veloped by the first author to create system quadru- 
ples of the discrete controller and discrete plant for 
FREQANAL.CMD. In general, the approach used is 
to combine feedback gains from a control-design pro- 
gram with other controller dynamics to form a state- 
space matrix that represents the controller transfer 
function. The plant is developed by combining indi- 
vidual continuous-system matrices (A,B,C,D) that 
have been previously constructed. These matrices 
are then assembled into a system quadruple and dis- 
cretized to form the discrete plant, which is used in 
most of the frequency-evaluation analyses. A second 
discrete plant, a combination of an airplane model 
with a Dryden wind-gust model, is also generated 
for use in the covariance analysis. All the various 
steps are accomplished by UDF’s that have been 
assembled to form SYSTEM.CMD. 

4.0. Time Analysis 

The time-analysis command files give the user a 
capability to evaluate control-system time responses 
as a function of signal commands or random exter- 
nal disturbances. The approach taken is similar to 
that described in section 3.0. Figure 9 contains the 
time-analysis flowchart. TIMEANAL.CMD executes 
other command files via menu selection to either pre- 
process data, build and simulate the control system 
to obtain time histories, or plot the generated data. 
The controller is assumed to be discrete, but the 
plant can be either discrete or continuous in the time 
simulations used with signal commands. An explicit 
solution is used to evaluate the all-discrete simula- 
tions, whereas a combination of a discrete solution 
and numerical integration is used for the hybrid time 
simulations. The time response to random external 


disturbances assumes a discrete plant that includes 
the disturbance model. 

BUILD is used to assemble the closed-loop 
control-system configurations, which consist of linear 
models of the plant and controller. An option exists 
to incorporate actuator saturation. The controller 
structure is general and allows for both feedback and 
feedforward compensation, which is also assembled 
in BUILD. 

The command files that can be selected via the 
time-analysis menu are as follows: 

1. TEVAL.CMD — computes time response to a 
signal-generator command 

2. TGUST.CMD — time response to a random 
input 

3. TPLOT.CMD — used for displaying time 
histories that are generated in both 
TEVAL.CMD and TGUST.CMD 

4. USERFILE — an optional file developed by the 
user 

It is possible for the general user to execute any of 
these files, as long as the appropriate variables are 
available. The' subsections that follow describe the 
contents of each command file and associated UDF 
in more detail. 

4.1. Time Evaluation 

All the time responses, except those due to a 
random external disturbance, are computed using 
TEVAL.CMD. This file serves as an executive for 
both discrete and hybrid closed- loop-system simula- 
tions. The hybrid closed-loop system has a discrete 
controller and a continuous plant. Menu selection 
allows the user to choose whether a limiter is to be 
used for actuator position saturation. If the choice is 
to use a limiter, then it is assumed that the variables 
for the upper and lower saturation bounds have been 
previously loaded into MATRIXx. 

TEVAL.CMD executes either DC1.CMD (if a dis- 
crete plant is to be used) or DC2.CMD (if a con- 
tinuous plant is to be used), as shown in figure 9. 
Both DCl.CMD and DC2.CMD build the closed- 
loop control system and connect the system to a sig- 
nal generator. After the total system is constructed, 
DC1.CMD or DC2.CMD runs the simulation. The 
output of these command files is a matrix with col- 
umn vectors that represent the output time histo- 
ries. DCl.CMD and DC2.CMD are nearly identi- 
cal. They both call GENC.FUN, which assembles 
the discrete generic controller, and they both exe- 
cute SIGNAL.CMD for the signal generator. The 
difference is that DCl.CMD calls DISCRETE.FUN, 
which assembles the discrete control system, whereas 
DC2.CMD calls HYBRID. FUN, which assembles the 
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hybrid control system. These files are each discussed 
in subsequent subsections. 

4*1*1 . Generic- controller structure . A 

generic-controller configuration is available for use 
within GENC.FUN. As shown within the dotted lines 
in figure 10, the structure allows the inclusion of 
a discrete feedback control C2 and a discrete feed- 
forward control, Cl or C3. BUILD is used to assem- 
ble and connect the blocks. Logic within GENC.FUN 
checks the quadruples SCI and SC3 to determine 
whether the corresponding Cl or C3 is to be con- 
nected. For example, if | |SC1| I 2 is zero, then Cl is 
not used and C3 is evaluated. If 1 1 SC3| (2 is also zero, 
then an error message is displayed, which indicates 
that there is no controller, and the program aborts. 
The feedback control quadruple that corresponds to 
C2 must be in the form 


SC2 = 


A B b B f 
C B b D fm 


( 8 ) 


where and T> b are the input matrices for the 
feedback measurements, and By and D y are the 
input matrices for the feedforward signals from C3. 
If Cl is used instead of C3, then matrices By and 
D f are not needed. If Cl is used as the feedforward 
control, then the output vector from Cl must be the 
same size as the output vector from C2, which is also 
the number of control signals going to the plant P. 

As shown in figure 11, various controller config- 
urations are possible with this structure. If Cl is a 
null matrix and C3 an identity matrix, then the con- 
figuration in figure 11(a) is obtained by making B^ 
and B y equal and making D b and D y equal. A neg- 
ative sign at the summing junction can be obtained 
by making By and Dy the negative of B b and D&, 
respectively. The configuration in figure 11(b) is ob- 
tained by making Cl a null matrix and letting C2 
have the general form of equation (8). Finally, the 
configuration shown in figure 11(c) is obtained by us- 
ing Cl and C2, with By and Dy not included in C2, 
and with C3 equal to zero. If C2 is a null matrix 
and Cl is an identity matrix, the command signal is 
input directly to the actuators, and the open-loop- 
plant time response can be obtained. 

4-1-2. Signal generator. SIGNAL.CMD is an 
executive command file for selecting a specific type 
of signal to be used for the simulation. A maximum 
of six signal generators of the same type can be 
assembled using BUILD. The simulation output time 
vector is also calculated in this file. First, the user 
is prompted for the maximum time TMAX of the 
simulation run. Next, a prompt appears for the user 
to decide whether the output time vector is to be 


limited to approximately 100 data points. If the 
response is “yes,” a second prompt then appears and 
requests the user to input the time step DELTIME 
between each element in the output time vector. The 
output time step T$ is then calculated to be some 
multiple of DELTIME by evaluating the equation 

T s = DELTIME x ROUND[(TMAX/DELTIME)/100] (9) 

where ROUND is a MATRIXx function used to 
compute the integer closest to the calculated value. 
The software does not allow this integer to be less 
than one. 

Equation (9) is used for both discrete and con- 
tinuous systems, but for a discrete system, the user 
should input the discrete sampling period (DELT) 
so that the output time step will be an even multiple 
of DELT. For a hybrid system, the user can input 
any number and Tg will be computed as shown in 
equation (9). If the user decides not to limit the out- 
put time vector to approximately 100 data points (an 
option available only for hybrid systems), the num- 
ber of data points is determined directly by dividing 
TMAX by DELTIME. A very large number of data 
points can be calculated if the user is not cautious. 

SIGNAL.CMD can handle multiple model cases. 
A menu option allows the user to select whether the 
same signal is to be used for all models or if different 
signals are to be used. The variable ANSW from this 
menu option is saved for use in the command files 
that build the specific signal generators. Another 
menu appears for the user to select the specific signal 
to be used. Two signal generators are currently 
available; one signal generator is for a pulse command 
(PULSE.CMD) and the second signal generator is for 
a ramp command (RAMP.CMD). Two of the signal 
sources available within BUILD, “pulse train” and 
“ramp,” are basic building blocks for PULSE.CMD 
and RAMP.CMD. 

PULSE.CMD can be used to create several differ- 
ent types of pulse signals, such as one pulse of vari- 
able width and amplitude, a step command signal, 
or a pulse train. A series of prompts allow the user 
to select the signal command paths for each run and 
to select the signal pulse width and the signal pulse 
amplitude. If different signals are to be used for the 
other models to be evaluated, these three prompts 
are repeated for each model. The start time is ini- 
tially defaulted to zero for all signals, but a menu 
option is available to change this time. Currently, all 
signals must be set to the same starting time, but 
this restriction may be removed in the future. 

There are three variables that the user inputs 
via prompts in RAMP.CMD. The user is required to 
select the signal-path number for each run, enter the 
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ramp slopes, and enter the ramp saturation value. 
All signal paths that are not used have amplitudes 
set to zero. The start time for the ramp command is 
always set to zero. If different signals are to be used 
for other models, the three prompts are repeated for 
each model. 

4-1-3. Control-loop construction . The re- 
maining task is to assemble the generic controller, 
plant, and limiter, if used, into a closed- loop con- 
trol system (SYS) and connect the signal gener- 
ator as shown in figures 12(a) and 12(b). The 
BUILD superblocks (a collection of blocks, ref. 8) 
created in GENC.FUN and SIGNAL.CMD are used 
in this assembly. The system is assembled in 
DISCRETE. FUN and HYBRID. FUN, with the 
superblock CGENC (fig. 12(b)) used in 

HYBRID. FUN and the superblock DGENC used in 
DISCRETE. FUN. DISCRETE. FUN is used for a 
discrete closed-loop control system, while 

HYBRID. FUN uses a continuous plant to form a hy- 
brid system. The simulation for the discrete system 
is valid only at the sampling intervals, but is more ef- 
ficient than the hybrid time simulation, since it takes 
less computer time. The hybrid system requires nu- 
merical integration for the plant equations of motion, 
but outputs can be observed between sampling times. 

Figure 13 shows a configuration for a hybrid 
simulation. The example includes a continuous plant 
with actuator saturation, a discrete controller, and a 
signal generator that can output either a pulse or a 
ramp with saturation. Outputs from both the plant 
and signal generator are used for display. 

One current restriction is that the vector of plant 
outputs must be composed of outputs from the ve- 
hicle, followed by outputs from the controls. The 
outputs for the controls can be either commands at 
the control inputs or control-surface positions if ac- 
tuator dynamics are included as part of the plant 
model. This restriction may be removed in a future 
modification. 

4.2. Time Response to External Disturbance 

The second option in the time analysis menu is 
TGUST.CMD (fig. 9), which evaluates the time re- 
sponse due to a random disturbance. Typically, the 
random disturbance considered by flight-control en- 
gineers is from wind gusts, but the computer code is 
general and allows any type of external-disturbance 
model. The user must supply a discrete plant sys- 
tem quadruple (SGDP) that is a combination of the 
original plant model and a disturbance model. There 
are certain input and output restrictions for SGDP 
that must be met. First, the number of inputs to 
the plant must equal the number of controller out- 
puts plus two additional inputs for the random dis- 


turbance, and the two disturbance inputs must be 
the last two input paths. Second, the last two out- 
puts from the plant must be from the disturbance 
model, since these outputs are later normalized to 
unity. If there is only one input disturbance of inter- 
est, the second disturbance path can always be used 
as a dummy input. The plant can have as many out- 
puts as desired. Connections from the plant to the 
controller are defined in the variable YVV. 

A typical example that a user might consider is a 
combined plant composed of a longitudinal airplane 
model and a Dryden wind-gust model with two out- 
puts, one output for the longitudinal axis and one 
output for the vertical axis of the airplane. Airplane 
time responses, including control-surface positions, 
can then be evaluated as a function of each separate 
input disturbance. Therefore, this computer code al- 
lows both axes to be evaluated during the same set 
of runs. 

As shown in figure 9, TGUST.CMD calls 
GUSTR.FUN. A menu allows selection of the in- 
put disturbance path, with only one input being 
selectable at a time for a given simulation run. 
GUSTR.FUN uses BUILD to construct the closed- 
loop control system. The only two inputs to this 
closed-loop system are the random disturbances. A 
white-noise disturbance of zero mean and normal 
distribution is input to the plant. A menu within 
GUSTR.FUN allows the user to select the default 
number of 1000 data points or to insert a differ- 
ent number of data points. Each data point is 
equivalent to one discrete sampling period, which 
is the conversion factor for the output time vector. 
The MATRIXx function DLSIM is used to gener- 
ate the simulation data. The standard deviation of 
the disturbance-model output (either of the last two 
outputs from the combined plant) is calculated first; 
all vectors in the output matrix are then normal- 
ized by this standard deviation. The output vectors 
are available for plotting, and the standard deviation 
and mean of all outputs are available in matrix form. 
(See appendix B.) As with the other analysis files, 
multiple models can be processed in the same run. 

4.3. Plotting 

The data generated in TEVAL.CMD and 
TGUST.CMD can be displayed by selecting 
TPLOT.CMD from the main time-analysis menu. 
There is one menu in TPLOT.CMD, and it offers 
a choice between two plot files named TEVAL.PLT 
and GUSTS. PLT. TEVAL.PLT is used to plot the 
time-history data generated in TEVAL.CMD, and 
GUST. PLT is used to plot the time-history data gen- 
erated in TGUST.CMD. 
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When TEVAL.PLT is selected, a display is shown 
with information such as the total number of data 
columns, number of runs, number of models used, 
number of plant outputs, number of controls, num- 
ber of command model outputs, and number of signal 
commands. Some of the displayed data are input to 
TEVAL.PLT, and some of the data are calculated 
within this command file. There is also a prompt for 
the user to enter a vector of variables desired for plot- 
ting. The user must know the order of data columns 
in terms of corresponding variables, and then must 
enter the column numbers to be plotted. However, 
the order of data columns is consistent with the order 
of outputs from the simulation model. The display is 
a helpful reminder as to the number of variables in- 
volved. The user can insert a display title by entering 
the desired characters within apostrophes. 

TGUST.PLT plots time-history data generated 
by TGUST.CMD. When this command file is se- 
lected, a display is shown with information similar 
to that described for TEVAL.PLT. The value for the 
disturbance-model output represents the two distur- 
bance paths that are used to generate the data. The 
user is prompted to enter a vector of variables to be 
plotted. Again, the user must know the order of the 
data columns in terms of the corresponding variables, 
but the display serves as a useful reminder. 

4.4. User File 

An optional user command file is available for 
the time analysis (fig. 9) that is similar to the user 
file described in section 3.3 for the frequency anal- 
ysis. When this option is selected, a prompt ap- 
pears for the user to enter a command-file name 
(within apostrophes). The user-file option has been 
included mainly as a convenience to allow the user a 
method of automatically processing data and gener- 
ating the required inputs for use in TEVAL.CMD and 
TGUST.CMD. This option does not have to be in- 
voked if the user desires to load the needed data into 
MATRIXx prior to executing TIMEANAL.CMD. A 
second use for this file could be to post-process data 
and save the data within a selected file. 

5.0. Integration With Control-Design Tools 

A goal of this integrated approach is to reduce the 
total overall design time by allowing the designer to 
quickly and accurately evaluate the control system 
by using several types of analysis. The idea is 
to combine separate control-design tools and these 
analysis tools in a manner that allows easy access 
to and between the design and analysis tools. The 
user-friendly control-system analysis tools described 
in this paper minimize the risk of human errors, since 


most of the analysis configurations are automatically 
assembled. 

An example of the integrated approach used by 
the authors is illustrated in figure 14. A user input 
data file is first loaded into MATRIXx for processing 
prior to executing a control-design program. Internal 
to the block labeled “PROCESS INPUT DATA” are 
multiple UDF’s that generate all data needed by the 
control-design program. As shown in figure 15, the 
UDF’s are combined to form a command file, and 
each UDF performs a specific task. The first task 
selects from a data base the plant model(s) to be used 
in the control design. All UDF’s have been developed 
to process multiple models automatically in the same 
run. The second task generates a reduced model by 
selecting the states and controls that are of interest. 
The third task selects the vehicle outputs that are to 
be used for control feedback. Other UDF’s have been 
developed that combine the vehicle and actuators to 
form a new plant, discretize the plant according to 
the selected sampling period, scale selected variables, 
convert state and control weightings to equivalent 
discrete weightings, and add both process and sensor 
noise variances. 

The control-design program presently being used 
is written in FORTRAN, but could be in any lan- 
guage. Since existing design programs are available, 
it seems reasonable that the newly developed tools 
interface with these programs. The data gener- 
ated by “PROCESS INPUT DATA” are saved on 
an ASCII file. These data can then be read by other 
programs. 

Data generated by a control-design program, typ- 
ically feedback gains (but could be any data such 
as eigenvectors and eigenvalues), are saved on a file 
which is then loaded back into MATRIXx. These 
data are then used by the analysis tools that have 
been discussed. The blocks shown within the dot- 
ted lines in figure 14 represent the external process 
performed by either the designer or some optimiza- 
tion process. After reviewing data generated by the 
analysis files, a decision is made to either accept the 
design or change some input design parameters. If 
the latter decision is made, the cycle is repeated. If 
the design is acceptable, a typical next step is to eval- 
uate the control law in a nonlinear simulation. 

6.0. Sample Problem 

This section contains descriptions and results of 
sample problems that are illustrated in detail in 
appendixes A and B. To make the sample problems 
easier to follow, the organization of appendixes A 
and B correspond to the organization of sections 6.1 
and 6.2, respectively. 
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6.1. Description of Frequency- Analysis 
Sample 

Appendix A contains the input data and the 
resulting screen output for the frequency analysis. 
Appendix B contains similar information for the time 
analysis. The sample plant is an airplane model 
with nine states (four vehicle states and five actuator 
states) and one control. 

6.1 .1. Input data. In appendix A, a discrete- 
plant quadruple (SDP) is shown with nine states 
(NS), one control, and three outputs, and a discrete- 
controller quadruple (SDC) is shown with two states 
(NSDC), three controls, and one output. Recall that 
the number of controller outputs must equal the num- 
ber of plant inputs, and the number of plant outputs 
must equal the number of controller inputs. The sam- 
pling period (DELT) is 0.031 sec (to be consistent 
with the corresponding real-time simulation model 
at Langley Research Center), and only one model is 
being evaluated. When more than one model is to be 
evaluated, SDP and SDC will be stacked by columns 
with the first model in the first set of columns. The 
quadruples (SDP and SDC) are used in the evalua- 
tions of singular- value robustness, singular- value loop 
transfer, Bode frequency response, and closed-loop 
eigenvalues. The steady-state covariance analysis re- 
quires a plant with an external-disturbance model. 
A third-order Dryden wind-gust model is combined 
with the plant (SDP) to evaluate the system response 
to external wind gusts in both the longitudinal and 
vertical directions. This new combined discrete-plant 
quadruple (SGDP) has 12 states (NSGDP), 3 con- 
trols, and 7 outputs. Two of the additional inputs 
and outputs are for the wind-gust model, and the 
other two additional outputs are for the actuator po- 
sitions. The controller and sampling period remain 
the same. 

6.1.2. Singular-value robustness. The 

singular- value multiplicative errors at the plant in- 
put and at the plant output are evaluated first. The 
frequency range selected is from 0.01 to 10 rad/sec 
with a density of 20 data points per decade. Since 
there is only one input, the structured and 
unstructured singular values are identical for this 
case. The structured singular value is selected for the 
outputs and new scale factors are calculated every 
fifth frequency point. The Davidon- Fletcher- Powell 
convergence algorithm is also selected. Outputs from 
the first two convergence passes are shown in appen- 
dix A as displayed on the computer terminal. The 
data point for the first convergence pass is calcu- 
lated twice to get good starting values. The data 
shown are the minimum singular value (SVALI), final 
step size (KAPPA), gradient of the two scale factors 


(GRADI), and the two scale factors (DOUT). The 
first path always has a scale factor of one; therefore, 
the number of scale factors is always one less than the 
number of paths. Figure 16 contains log-log plots of 
the minimum singular value at the plant input (SIGI) 
and the minimum singular value at the plant output 
(SIGO). 

6.1.3. Singular-value loop transfer and 
Bode frequency response. Both the singular- 
value loop transfer and Bode frequency response 
are calculated over a revised frequency range of 
0.1 to 100 rad/sec with a density of 20 points per 
decade. Figure 17 contains log-log plots of the max- 
imum singular value at the plant input (SVLTI) 
and the maximum singular value at the plant out- 
put (SVLTO). There are two plots that illustrate 
the Bode frequency response. Figure 18 contains 
semilog frequency- response plots of gain in decibels 
and phase angle. A second illustration is shown in 
figure 19, which is a Nyquist polar plot of ampli- 
tude versus phase angle. The amplitude scale for 
the Nyquist plot is dimensionless, as opposed to dB, 
since the radius of the MATRIXx polar plot must be 
positive. 

6.1.4 • Steady-state covariance. Typical 
results from the steady-state covariance analysis are 
illustrated in the bar chart in figure 20. Selections 1 
(XURMS) and 2 (YURMS) are for the state and 
output RMS values, respectively, with a normalized 
wind-gust disturbance in the first disturbance path. 
States 4, 5, and 6 are plotted in the top portion of 
figure 20, and outputs 3, 4, and 5 are plotted in the 
bottom portion of the figure. 

6.1.5. Closed-loop eigenvalues. The final 
evaluation is the closed- loop eigenvalue test. The 
tabular eigenvalue data (EVALCL) are calculated 
and remain on file for observation after exiting the 
FREQANAL.CMD menu but prior to exiting 
MATRIXx. The table is set up with the complex 
discrete eigenvalues in the first two columns, the 
complex continuous eigenvalues in the second two 
columns, and the damping ratios in the real part of 
the last two columns. 

6.1.6. Additional output data. Other out- 
put data (in addition to EVALCL) include the phase 
margin (PHASEM) and gain margin (GAINM). The 
four columns of PHASEM represent the phase mar- 
gin in degrees, frequency in rad/sec, model num- 
ber, and control loop number. The first column 
of GAINM is the gain margin in dB; the remain- 
ing three columns represent the same quantities as 
in PHASEM. The RMS data (XURMS, YURMS, 
XWRMS, YWRMS) for the two external-disturbance 
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paths are also shown. XURMS and XWRMS have 
14 elements corresponding to the 12 states in SGDP 
plus 2 states in SDC; YURMS and YWRMS have 
7 elements corresponding to the 7 outputs from 
SGDP. 

The example described is for one model. When 
multiple models are evaluated, the output vectors are 
stacked by columns; the first set of vectors represents 
the first model. As an example, for the case of q 
models and p outputs, the first p columns of the 
output matrix will be from the first model, the second 
p columns will be from the second model, and the 
total number of columns in the output matrix will 
be the product of p and q. This same data stacking 
approach is used in the time analysis. 

6.2. Description of Time-Analysis Sample 

The input data and resulting screen output for 
the time-analysis sample described in this subsection 
are shown in appendix B. The sample plant is an 
airplane model with 9 states (4 vehicle states and 
5 actuator states) and 2 controls. 

6-2-1. Input data. Input data for the time 
analysis shown in appendix B are slightly more com- 
plicated than those needed for the frequency analysis. 
The generic controller (SCI, SC2, SC3) has dynamic 
feedback (SC2) with two states (NSC2), two out- 
puts, and six inputs. The six inputs are partitioned 
into three inputs for the feedback signals from the 
plant and three identical inputs of opposite sign for 
the commands from the signal generator (columns 6 
to 8 of SC2 are the negative of columns 3 to 5). This 
set of inputs is equivalent to a summing junction be- 
tween the feedback signals and the command signals, 
as shown in figure 11(a). The feedforward control is 
just a unity gain matrix (SC3). Note that SCI and 
NSC1 must be set to zero if they are not used. 

The plant selected for the sample problem is con- 
tinuous (SCP) with nine states (NSCP), two con- 
trols, and seven outputs. The seven outputs include 
five outputs from the vehicle and two outputs for the 
control positions. Three of the vehicle outputs are 
fed back to SC2, and the other outputs are used only 
for display. Three quantities that must be defined 
are located in variable DIMS. DIMS (9) contains 
the sampling period of 0.031 sec, DIMS(10) contains 
the number of vehicle outputs, and DIMS(ll) is the 
starting point for a row vector that defines the three 
vehicle outputs that are fed back. Three other addi- 
tional inputs required are the number of models used 
(NMOD) and row vectors that contain the lower and 
upper saturation values of the actuators (LBOUND, 
UBOUND). These last two variables are required; if 
a limiter is not being used, the variables can be set 
to the scalar zero. 


The input data for the time response to exter- 
nal disturbances (TGUST.CMD) consist of a dis- 
crete quadruple for the combined plant (SGDP) with 
12 states (NSGDP), 4 inputs, and 9 outputs. The in- 
puts to SGDP consist of two controls and two random 
white-noise signals with normal distribution. The 
discrete controller (SDC) has two states (NSDC), two 
outputs that go to the plant, and three measurement 
inputs from the plant (YVV). 

6-2.2. Time evaluation. In the first example, 
TEVAL.CMD is selected. A continuous plant with 
actuator saturation is being used for the 10-sec sim- 
ulation with approximately 100 data points of out- 
put at time intervals which are a multiple of DELT. 
A ramp signal generator with a slope of 0.2 and a 
saturation level of 1 is the command to the first in- 
put path; the other two nonselected input paths have 
default input commands of zero. 

After this simulation run is completed, a second 
run is selected for the same maximum time and num- 
ber of data points. The command to the first path is 
a step with an amplitude of 100; implementation of 
this step is accomplished by selecting a pulse and set- 
ting the pulse width equal to the maximum time of 
the simulation run (10 sec). Defaults include setting 
the start time of the pulse to zero, setting the am- 
plitude of the nonselected signal generators to zero, 
and setting the pulse frequency to a small number. 

An option exists to save or delete the files within 
BUILD. It is generally advisable to delete these files, 
since they will remain with other data that the user 
may desire to save. If the BUILD files are kept, it is 
then possible to edit these files, after finishing with 
TIMEANAL.CMD, and to run special cases. 

The plotting option is selected with the option 
of a hard-copy output on file MATLASER.DAT. 
A split screen is selected with four variables for 
the upper plot and three variables for the lower 
plot. Prior to the selection of the variable, a display 
shows the user that there are 22 columns of data 
for the two runs and that the data for each run 
are distributed as 5 columns for the plant (airplane) 
outputs, 2 columns for the control outputs from the 
plant, 3 columns for the model outputs, and 1 column 
for the signal-generator output. The simulation title 
REGULATOR is attached to the plot as shown in 
figure 21. Variable names cannot be printed within 
a key, but the line style of the curves always appears 
in a known order. 

6.2.3 . Time response to external distur- 
bance. TGUST.CMD is selected as another exam- 
ple. First, the random disturbance for the first input 
path (channel 1) is selected with the default value 
of 1000 data points. After this run is complete, the 
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second input path (channel 2) is selected. The ex- 
ample shows that the number of data points can be 
changed for each time-response calculation. An ex- 
ample containing upper and lower plots is shown in 
figure 22. Two of the nine available variables are plot- 
ted in the upper part of the figure, and one variable 
is shown in the lower part of the figure. 

6.2-4- Additional output data . After quitting 
TIMEANAL.CMD, the generated data are available 
in memory. Data that are needed must be saved prior 
to exiting MATRIXx, or they will be lost. Examples 
of statistical data for the external-disturbance input 
paths (channels 1 and 2) are illustrated by GUSTAT 
and GWSTAT, respectively. The first row of each 
variable represents the standard deviation of the 
time-history data for each displayed output (defined 
by columns), and the second row represents the 
mean of the time-history data. All data are for a 
normalized input disturbance of 1. 

7.0. Concluding Remarks 

One objective of this work is to enable the user 
to assemble and analyze various control-system con- 
figurations within the MATRIXx software package. 
Implementation of a high-level system approach to 
control- law analysis has been accomplished by com- 
bining basic MATRIXx functions with user-defined 
functions and command files. Programming features 
within MATRIXx and its BUILD module allow var- 
ious control-system analysis configurations to be as- 
sembled automatically after the user selects the anal- 
ysis to be made. Embedded within these tools are 
menu-driven options that require minimal input data 
and interaction by the user. This approach minimizes 
the mistakes by the designer and the time required 
to perform the various analyses. Key attributes of 
these control analysis tools include simplicity, user- 
friendliness, flexibility, and the ability to trade off 
between time and accuracy. 

Frequency-response analysis tools and time- 
response analysis tools are both available. All these 
analysis routines have been developed to handle 
multiple models automatically in one computer run. 
This feature allows, in minimum time, compari- 
son of different control systems at different flight 
conditions. 

The only inputs required for the frequency evalu- 
ations are discrete state-space quadruples for a plant 
and a controller, the corresponding number of states 
for each quadruple, and the discrete sampling pe- 
riod. Using these inputs, five different evaluations 
can be performed. The input data are minimized for 
simplicity. 


A vector of frequency data points, based upon 
user inputs which define the frequency range of inter- 
est and the density of data points within that range, 
is automatically created by the program. This fre- 
quency vector can be easily changed for each analy- 
sis. The frequency vector can be composed of sev- 
eral subvectors, each with a different density of data 
points. This feature is important when evaluating 
a control system with a lightly damped eigenvalue, 
because additional data points can be calculated in 
the frequency range where the control-response vari- 
able changes rapidly. In contrast, fewer data points 
are needed in frequency ranges with slowly varying 
control responses. This flexibility results in reduced 
time to perform the analysis, but still allows the im- 
portant response features to be accurately analyzed. 

The frequency evaluations are applicable to both 
single input-output systems and multiple input- 
output systems. One frequency-evaluation test 
developed in recent years for multiple input-output 
systems is the structured singular value for a 
multiplicative perturbation at either the plant input 
or the plant output. Theoretically, a new scaling ma- 
trix must be computed at each new frequency. Since 
the scaling matrix must be solved iteratively, the pro- 
cess is time-consuming. The user has the option of 
reducing the number of frequencies at which the scal- 
ing matrix is calculated; this option decreases the 
computation time with the penalty of a slight de- 
crease in accuracy. It is also possible to compute the 
unstructured singular- value frequency response. 

A generic controller has been developed for use 
in simulations that use signal-generator commands. 
The structure allows the incorporation of a dynamic 
feedback controller and two options for a dynamic 
feedforward controller. Several different control- 
system configurations can be assembled and evalu- 
ated, including configurations that allow signal com- 
mands to be input directly to the actuators. 

The generic controller must be in discrete form, 
but the plant can be either discrete (for an all- 
discrete simulation) or continuous (for a hybrid sim- 
ulation). The all-discrete simulation is valid only 
at the sampling intervals, but it is more efficient 
than the hybrid simulation because it takes less com- 
puter time. The hybrid system requires numerical 
integration for the plant equations of motion, but 
outputs can be observed between sampling times. 
With either type of simulation, the user may choose 
whether a limiter is to be used for actuator position 
saturation. 

The time-analysis command files give the user 
the capability to evaluate control-system time re- 
sponses as a function of signal commands or ran- 
dom external disturbances. Presently, either ramp 
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or pulse signal generators are available for command 
inputs, but others may be added in the future. The 
pulse command can be either a single pulse, a step, 
or a pulse train, and the ramp may have satura- 
tion included. Time responses and statistical data 
(mean and standard deviation) are generated when 
an external-disturbance input is used. The external 
disturbance is white noise with zero mean and nor- 
mal distribution. 

Options exist for the inclusion of user- developed 
command files. The user-file selection has been 
included mainly as a convenience to allow the user 
a method of automatically preprocessing data. This 
option does not have to be invoked if the user has 
already loaded the required data into MATRIXx 
prior to executing the selected command file. 

Plot command files are available for both fre- 
quency and time evaluations. Singular-value 
frequency-response data are displayed on log-log 
plots, Bode frequency-response data are displayed 
on semilog plots, and steady-state root-mean-square 
data due to an external disturbance are displayed in 
bar charts. Time-history outputs are displayed on 
linear scaled plots. 

A plan has been demonstrated for integration of 
the analysis tools with control-design tools; this plan 
allows easy access to and between the two processes. 
This integrated approach allows the designer to easily 
change design parameters, run the control-synthesis 
program, and then perform various analyses. Total 
time for this complete cycle is significantly reduced 
relative to an approach which has few integrated 
tools. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
January 11, 1989 
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APPENDIX A 


DATA FOR FREQUENCY-ANALYSIS SAMPLE 

The "Sample Problem" section (6.0) describes the input data needed for 
the frequency analysis. These data and the output as displayed on the screen 
are shown below. For clarity, the titles, section numbers, and user responses 
to prompts are in boldface print. Also, explanations are located on the right side 
of the page. 


A.1 INPUT DATA 


SDP 


COLUMNS 1 
9.9858D-01 
-2.2064D-05 
-3.0376D-06 
-4.8550 D-08 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
6.6066D-03 


THRU 6 
-4.9898 D-01 
9.9029D-01 
-2.4344D-02 
-3.7840D-04 
0.0000D+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
1.0000D+00 
O.OOOOD+OO 
2.8002D+00 


COLUMNS 7 THRU 10 


-1.9398D-02 
-3.5455D-04 
-5.5898D-03 
-1.6246D-04 
0.0000D+00 
O.OOOOD+OO 
-1 .2748 D-01 
-7.0635D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
2.651 5D-01 


-2.3460D-02 
-3.4715D-04 
-6.7448D-03 
-1.1376D-04 
O.OOOOD+OO 
O.OOOOD+OO 
7.0635D-02 
-1 .2748 D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


NS 

9. 

SDC = 

8.0880D-01 

O.OOOOD+OO 

1.0000D+00 


-7.8095D-01 

1.0000D+00 

O.OOOOD+OO 


-2.3061 D-02 
3.0419D-02 
9.9263 D-01 
3.0888D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
1.0000D+00 
9.7528D-02 


-1.1230D-02 

-2.1587D-04 

-5.1480D-03 

-9.2150D-05 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

3.9455D-01 

O.OOOOD+OO 

O.OOOOD+OO 

5.7684D-02 


-1.0667D+00 

3.1000D-02 

O.OOOOD+OO 


-9.9636D-01 
1.1012D-05 
1.5615D-06 
1.0000D+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
2.8682D-1 3 


4.1398D-03 
5.3172D-05 
1.2374D-03 
1.1557D-05 
-6.3848D+00 
-1.4840D+01 
-1.0190D-01 
-1.2554D-01 
-2.8471 D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-2.1447D+00 

3.1000D-02 

O.OOOOD+OO 


-6.8684D-04 
-8.9517D-06 
-1.9730D-04 
-2.1021 D-06 
3.2457D-01 
-5.3951 D-01 
1.1730D-02 
5.8018D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


1.2979D-01 

O.OOOOD+OO 

O.OOOOD+OO 


-1.9073D-04 
-2.3042D-06 
-5.4769D-05 
-3.9946D-07 
5.3951 D-01 
3.2457D-01 
6.3442D-03 
8.0366D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


16 


NSDC 


2 . 

DELT = 


3.1000D— 02 


SGDP = 


COLUMNS 1 
9.9858D-01 
-2.2064D-05 
-3.0376D-06 
-4.8550D-08 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
6.6066D-03 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 


IHRU 6 
-4.9898D-01 
9.9029D-01 
-2.4344D-02 
-3.7840D-04 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
1.0000D+00 
0.0000D+00 
2.8002D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 


-2.3061 D-02 
3.0419D-02 
9.9263D-01 
3.0888D-02 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
1.0000D+00 
9.7528D-02 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 


-9.9636D-01 
1.1012D-05 
1.5615D-06 
1.0000D+00 
0.0000D+00 
0.0000D+00 
O.OOOOD+OO 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
2.8682D-1 3 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 


-6.8684D-04 
-8.9517D-06 
-1 .9730D-04 
-2.1 021 D-06 
3.2457D-01 
-5.3951 D-01 
1.1730D-02 
5.8018D-03 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
O.OOOOD+OO 
0.0000D+00 
0.0000D+00 
O.OOOOD+OO 
0.0000D+00 
O.OOOOD+OO 
O.OOOOD+OO 


COLUMNS 7 
-1.9398D-02 
-3.5455D-04 
-5.5898D-03 
-1.6246D-04 
O.OOOOD+OO 
O.OOOOD+OO 
— 1 . 2748 D— 01 
-7. 0635 D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
2.6515D-01 
2.5144D+01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


rmu 12 

-2. 3460 D-02 
-3.4715D-04 
-6.7448D-03 
-1.1376D-04 
O.OOOOD+OO 
O.OOOOD+OO 
7.0635D-02 
-1 .2748 D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-1.1230D-02 

-2.1587D-04 

-5.1480D-03 

-9.2150D-05 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

3.9455D-01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

5.7684D-02 

O.OOOOD+OO 

3.0000D+01 

O.OOOOD+OO 

O.OOOOD+OO 


1 .4204D-03 
2.2006D-05 
3.0294D-06 
4.8464D-08 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
9.9475D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
1.0000D+00 
O.OOOOD+OO 


1.6756D-03 
3.2595D-05 
8.1752D-05 
1 .271 8D-06 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
9.9475D-01 
-5.2335D-05 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
3.3670D-03 


-1.9073D-04 
-2.3042D-06 
-5.4769D-05 
-3.9946D-07 
5.3951 D-01 
3.2457D-01 
6.3442D-03 
8.0366D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


4.41 28D-08 
8.4718D-1 0 
2.1547D-09 
2.231 6D-1 1 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
5.2335D-05 
9.9475D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+O'' 
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COLUMNS 13 
4.1398D-03 
5.31 72D-05 
1.2374D-03 
1.1557D-05 
-6.3848 D+00 
-1.4840D+01 
-1.0190D-01 
-1.2554D-01 
-2.8471 D-02 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 

NSGDP = 

12 . 


THRU 15 
1.2858D-05 
1.9912D-07 
2.8235D-08 
2.9631 D-10 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
1.8013D-02 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 


1.8572D-05 
3.5655D-07 
9.0683D-07 
9.3896D-09 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
2.2037D-02 
-9.3255 D-01 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 


A. 2 SINGULAR-VALUE ROBUSTNESS 


<> EXEC ('MUGDIR:FREQANAL.CMD') 
SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXECUSERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 1 
ANALYSIS: 

1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 

SELECT> 1 
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OMEG MENU: 

1 CHANGE OMEG 

2 REVIEW OMEG 

3 QUIT 

SELECT> 1 

NEW OMEG FREQ RANGES<MIN,MAX1,‘**,MAXN> [.01 10] 
N 


1 . 

PTS/DECADE FOR EACH RANGE(NOM=2u)<1,2,***,N> [20] 
OMEG MENU: 

1 CHANGE OMEG 

2 REVIEW OMEG 

3 QUIT 

SELECT> 3 
ANALYSIS: 

1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 

SELECT> 2 


OPTION 2. (SINGULAR VALUES - MULTIPLICATIVE ERROR) 
<OMEGA,SIG,SVAL,DSCALED,GRA>=SSVIO(S,NS,OMEG,DELT) 
<SIMB,NSIMB,SOM,NSOM>=SV(SDP,NS,SDC,NSDC,DELT) 

PLEASE CHOOSE: 

1 INPUTS ONLY 

2 OUTPUTS ONLY 

3 BOTH 

SELECT> 3 

NCOUNT SELECTION NCOUNT IS THE ITERATION RATE FOR 

CALCULATING NEW SCALE FACTORS. 

1 KEEP DEFAULT (5) 

2 CHANGE NCOUNT 

3 SELECT NCOUNT=0 


SELECT> 1 


FUNCTION SELECTION 


1 DEFAULT (SVDFP) 

2 SELECT SVGNR 

SELECT>1 

COMPLETED EXAMPLE FOR 1 INPUT 

NCOUNT SELECTION 

1 KEEP DEFAULT (5) 

2 CHANGE NCOUNT 

3 SELECT NCOUNT=0 

SELECT> 1 

FUNCTION SELECTION 

1 DEFAULT (SVDFP) 

2 SELECT SVGNR 

SELECT> 1 
SVALI = 

1 .0090 
KAPPA = 

1. 

GRADI = 

- 0.0012 

-0.0157 

DOUT = 

9.0403 

0. 0896 

COMPLETED FIRST POINT, INITIAL CALCULATION 

SVALI = 

1 .0089 
KAPPA = 

1. 

GRADI = 

-0.0013 

-0.0004 
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DOUT 


9.0614 

0.1473 


I 


1 . 

COMPLETED FIRST POINT, FINAL CALCULATION 

SVALI = 

1.0166 
KAPPA = 

1 . 

GRADI = 

1.0D-03* 

-0.8142 

-0.6040 

DOUT = 

9.0721 

0.1379 


2 . 


COMPLETED SECOND POINT 


COMPLETED EXAMPLE FOR 3 OUTPUTS 


SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXECUSERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 2 
AMP =AMPL(DB) 


PICTURE TAKING: 


1 EVERYTHING 

2 SOMETHINGS 

3 NOTHING 

SELECT> 2 


Available devices are: 

1 VT220 

2 VT240/Monochrome 

3 Tektronix 4014 

4 MATLN03 

SELECT> 4 

ENTER GRAPHICS FILENAME [MATLASER.DAT]: CARRIAGE RETURN 
PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 1 


PLOT MENU 

1 WHOLE 

2 1/2 UP/LOW 

3 1/2 UP. 1/4 LOW. 

4 1/4 UP. 1/2 LOW. 

5 1/4 QUAD 

6 BODE DEFAULT 

7 QUIT 

SELECT> 2 
MEASURED OUTPUTS 
NY 
3. 

CONTROL COMMANDS 
NBU 
1 . 
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SELECT FOR UPPER 


CHOOSE: 

1 SIGI 

2 SIGO 

3 QUIT 


SELECT> 1 


ENTER VECTOR OF MODELS DESIRED [1] 
SELECT FOR LOWER 


TOP HALF PLOTTED (FIG. 16) 


CHOOSE: 


1 SIGI 

2 SIGO 

3 QUIT 

SELECT> 2 

ENTER VECTOR OF MODELS DESIRED [1] 

LOWER HALF PLOTTED (FIG. 16) 


TAKE PICTURE NOW: 

1 YES 

2 NO 

SELECT> 1 
PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 5 


A.3 SINGULAR-VALUE LOOP TRANSFER AND BODE FREQUENCY RESPONSE 


SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXEC USERFILE 

4 EXEC AJO USERFILE 

5 QUIT 


SELECT> 1 



ANALYSIS: 


1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 

SELECT> 1 
OMEG MENU: 

1 CHANGE OMEG 

2 REVIEW OMEG 

3 QUIT 

SELECT> 1 

NEW OMEG FREQ RANGES<MIN,MAX1,***,MAXN> [.1 100] 
N 


1 . 

PTS/DECADE FOR EACH RANGE(NOM=20)<1,2,***,N> [20] 
OMEG MENU: 

1 CHANGE OMEG 

2 REVIEW OMEG 

3 QUIT 

SELECT> 3 


A.3.1 SINGULAR-VALUE LOOP TRANSFER 


ANALYSIS: 

1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 

SELECT> 3 

COMPLETED SINGULAR VALUE LOOP TRANSFER CALCULATION 
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SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXECUSERFILE 

4 EXECAJOUSERFILE 

5 QUIT 

SELECT> 2 
AMP =AMPL(DB) 

PICTURE TAKING: 

1 EVERYTHING 

2 SOMETHINGS 

3 NOTHING 

SELECT> 2 
PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 2 
PLOT MENU 

1 WHOLE 

2 1/2 UP/LOW 

3 1/2 UP. 1/4 LOW. 

4 1/4 UP. 1/2 LOW. 

5 1/4 QUAD 

6 BODE DEFAULT 

7 QUIT 

SELECT> 2 
MEASURED OUTPUTS 
NY = 

3. 

CONTROL COMMANDS 
NBU = 


1 . 


SELECT FOR UPPER 


CHOOSE: 

1 SVLTI 

2 SVLTO 

3 QUIT 

SELECT> 1 


ENTER VECTOR OF MODELS DESIRED <ON A SINGLE PLOT> [1] 

ENTER SV(S) TO BE PLOTTED [1] 

TOP HALF PLOTTED (FIG. 17) 


SELECT FOR LOWER 
CHOOSE: 

1 SVLTI 

2 SVLTO 

3 QUIT 


SELECT> 2 

ENTER VECTOR OF MODELS DESIRED <ON A SINGLE PLOT> [1] 

ENTER SV(S) TO BE PLOTTED [1] 

LOWER HALF PLOTTED (FIG. 17) 


TAKE PICTURE NOW: 

1 YES 

2 NO 

SELECT> 1 
PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 5 
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A.3.2 BODE FREQUENCY RESPONSE 


SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXEC USERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 1 
ANALYSIS: 

1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 


COMPLETED BODE FREQUENCY RESPONSE CALCULATION 
SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXEC USERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 2 
AMP =AMPL(DB) 

PICTURE TAKING: 

1 EVERYTHING 

2 SOMETHINGS 

3 NOTHING 

SELECT> 2 
PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 


SELECT> 3 


PLOT MENU 


1 WHOLE 

2 1/2 UP/LOW 

3 1/2 UP. 1/4 LOW. 

4 1/4 UP. 1/2 LOW. 

5 1/4 QUAD 

6 BODE DEFAULT 

7 QUIT 

SELECT> 2 

CONTROL COMMANDS 
NBU = 

1 . 

SELECT FOR UPPER 
CHOOSE: 

1 DB 

2 PHASE 

3 NYQUIST 

4 QUIT 

SELECT> 1 

ENTER VECTOR OF MODELS DESIRED <ON A SINGLE PLOT> [1] 

ENTER THE CHANNEL TO BE GRAPHED [1] 

TOP HALF PLOTTED (FIG. 18) 


SELECT FOR LOWER 
CHOOSE: 

1 DB 

2 PHASE 

3 NYQUIST 

4 QUIT 

SELECT> 2 

ENTER VECTOR OF MODELS DESIRED <ON A SINGLE PLOT> [1] 

ENTER THE CHANNEL TO BE GRAPHED [1] 

LOWER HALF PLOTTED (FIG. 18) 


TAKE PICTURE NOW: 

1 YES 

2 NO 

SELECT> 1 
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PLOT SELECT: 


1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 3 
PLOT MENU 

1 WHOLE 

2 1/2 UP/LOW 

3 1/2 UP. 1/4 LOW. 

4 1/4 UP. 1/2 LOW. 

5 1/4 QUAD 

6 BODE DEFAULT 

7 QUIT 

SELECT> 1 

CONTROL COMMANDS 
NBU = 

1 . 

FOR SINGLE PLOT 
CHOOSE: 

1 DB 

2 PHASE 

3 NYQUIST 

4 QUIT 

SELECT> 3 

ENTER VECTOR OF MODELS DESIRED <ON A SINGLE PLOT> [1] 

ENTER THE CHANNEL TO BE GRAPHED [1] 

CALULATING THE AMPLITUDE AT THIS TIME. 

NYQUIST POLAR PLOT COMPLETED (FIG. 19) 


TAKE PICTURE NOW: 

1 YES 

2 NO 


SELECT> 1 


PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 5 


A.4 STEADY-STATE COVARIANCE 


SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXECUSERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 1 
ANALYSIS: 

1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 

SELECT> 5 


COMPLETED STEADY-STATE COVARIANCE CALCULATION 


SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXEC USERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 2 
AMP =AMPL(DB) 

PICTURE TAKING: 

1 EVERYTHING 

2 SOMETHINGS 

3 NOTHING 

SELECT> 2 
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PLOT SELECT: 


1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 4 
PLOT MENU 

1 WHOLE 

2 1/2 UP/LOW 

3 1/2 UP. 1/4 LOW. 

4 1/4 UP. 1/2 LOW. 

5 1/4 QUAD 

6 BODE DEFAULT 

7 QUIT 

SELECT> 2 
MEASURED OUTPUTS 
NY 
3. 

GUSTS INPUTS = 2 
SELECT FOR UPPER 
CHOOSE: 

1 XURMS 

2 YURMS 

3 XWRMS 

4 YWRMS 

5 QUIT 

SELECT> 1 

ENTER VECTOR OF MODELS DESIRED [1] 

SELECT VECTOR OF ELEMENTS DESIRED [4 5 6] 

TOP HALF PLOTTED (FIG. 20) 


SELECT FOR LOWER 
CHOOSE: 

1 XURMS 

2 YURMS 

3 XWRMS 

4 YWRMS 

5 QUIT 


SELECT> 2 


ENTER VECTOR OF MODELS DESIRED [1] 

SELECT VECTOR OF ELEMENTS DESIRED [3 4 5] 

LOWER HALF PLOTTED (FIG. 20) 


TAKE PICTURE NOW: 

1 YES 

2 NO 

SELECT> 1 
PLOT SELECT: 

1 SV 

2 SVLT 

3 BFR 

4 GUST 

5 QUIT 

SELECT> 5 


A.5 CLOSED-LOOP EIGENVALUES 


SELECT COMMAND FILE 

1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXEC USERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 1 
ANALYSIS: 

1 CHANGE OMEG 

2 S.V. MULT ERR 

3 S.V. LOOP TRAN. 

4 BODE FREQ RESP 

5 SSCOV 

6 CLOSED LOOP EIG 

7 QUIT 

SELECT> 6 

COMPLETED CLOSED-LOOP EIGENVALUES CALCULATION 
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SELECT COMMAND FILE 


1 EXEC FREQEVAL.CMD 

2 EXEC PLOT.CMD 

3 EXEC USERFILE 

4 EXEC AJO USERFILE 

5 QUIT 

SELECT> 5 


A.6 ADDITIONAL OUTPUT DATA 


EVALCL = 


-1.2610D-01 

-1.2610D-01 

3.9387D-01 

3.2015D-01 

3.2015D-01 

9.1988D-01 

9.1988D-01 

9.7906D-01 

9.7906D-01 

9.9960D-01 

9.9960D-01 

PHASEM = 


-7.0068D-02i 
+7.0068D-02i 
—8.181 7D— 1 8i 
+5.301 0D-01I 
-5.301 OD-Oli 
-9.4606D-02i 
+9.4606D-02I 
-1.0570D-02i 
+1.0570D-02i 
-4.0888D-03i 
+4.0888D-03i 


-6.2456D+01 
-6.2456D+01 
-3.0056D+01 
-1.5458D+01 
-1.5458D+01 
-2.5242D+00 
-2.5242D+00 
-6.8090D-01 
-6.8090D-01 
-1 .2555 D-02 
-1.2555D-02 


-8.4981 D+Oli 
+8.4981 D+Oli 
-6.7009D-16i 
+3.3145D+01 i 
-3.3145 D+Oli 
-3.3060D+00i 
+3.3060D+00i 
-3.4825D-0H 
+3.4825D-0U 
-1.3195D-0U 
+1.3195D-0U 


5.9220D-01 
5.9220D-01 
1.0000D+00 
4.2268D-01 
4.2268D-01 
6.0686D-01 
6.0686D-01 
8.9031 D-01 
8.9031 D-01 
9.4722D-02 
9.4722D-02 


5.0008D+01 
GAINM = 

-1 .6040D+01 


3.0826D+00 

1.1110D+01 


1.0000D+00 

1.0000D+00 


1.0000D+00 

1.0000D+00 


XURMS = 

1.6452D+00 
9.9085 D-04 
9.1346D-04 
6.2951 D-03 
2.6258D-02 
1.4032D-02 
5.1980D-04 
3.7986D-04 
7.4870D-05 
1.0000D+00 
1.2548D-15 
2.1641D-15 
1.5925D-03 
9.6856D-04 


+0.0000D+00i 

+0.0000D+00i 

+O.OOOOD+OOi 

+0.0000D+00I 

+0.0000D+00i 

+0.0000D+00i 

+0.0000D+00i 

+0.0000D+00i 

+0.0000D+00i 

+0.0000D+00i 

+0.0000D+00i 


YURMS 


9.9085D-04 
9.1346D-04 
9.0288 D-03 
1.3070D-02 
2.2461 D-03 
1.0000D+00 
4.2248 D-1 8 

XWRMS = 

5.9138D+02 
4.8425D-01 
4.0372D-01 
2.4170D+00 
4.0363 D+01 
2.1 609D+01 
7.9894D-01 
5.8395D-01 
1.1496D-01 
6.5833D-12 
2.9700D+02 
1.5378D+04 
2.4466D+00 
8.9393D-01 

YWRMS = 

4.8425D-01 

4.0372D-01 

3.2727D+00 

2.0088D+01 

3.4487D+00 

6.5833D-12 

1.0000D+00 


APPENDIX B 


DATA FOR TIME-ANALYSIS SAMPLE 


The "Sample Problem" section (6.0) describes the input data needed for the 
time analysis. These data and the output as displayed on the screen are shown 
below. For clarity, the titles, section numbers, and user responses to prompts are 
in boldface print. Also, explanations are located on the right side of the page. 


B.1 INPUT DATA 


SCI = 

0 . 

NSC1 = 

0 . 

SC2 = 

COLUMNS 1 THRU 6 

8.0880D-01 -7.8095 D-01 -1.0667D+00 -2.1447D+00 1.2979D-01 1.0667D+00 

0.0000D+00 1.0000D+00 3.1000D-02 3.1000D-02 0.0000D+00 -3.1000D-02 

-8.2059D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 

-1.4107D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 

COLUMNS 7 THRU 8 
2.1447D+00 -1.2979D-01 
-3.1000D-02 0.0000D+00 

0.0000D+00 0.0000D+00 

0. 0000D+00 0.0000D+00 

NSC2 = 

2 . 

SC3 = 

1 . 0 . 0 . 

0 . 1 . 0 . 

0 . 0 . 1 . 

NSC3 = 

0 . 
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SC P 


COLUMNS 1 THRU 6 

-4.6153D-02 -1.6189D+01 2.8067D-13 

-7.1404D-04 -3.0264D-01 9.8946D-01 

-1.0716D-04 -7.9189D-01 -2.2643D-01 
0.0000D+00 0.0000D+00 1.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

1.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 1.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 1.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

6.6066D-03 2.8002D+00 9.7528D-02 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

COLUMNS 7 THRU 11 

-2.8605D+00 0.0000D+00 -5.6114D-01 

-2.8658D-02 0.0000D+00 -6.2345D-03 

-8.21 24D-01 0.0000D+00 -2.5625D-01 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

-6.2127D+01 8.5020D+01 0.0000D+00 

-8.5020D+01 -6.21 27D+01 0.0000D+00 

0.0000D+00 0.0000D+00 -3.0000D+01 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 O.OOOOD+OO 

0.0000D+00 0.0000D+00 0.0000D+00 

2.6515D-01 0.0000D+00 5.7684D-02 

2.5144D+01 0.0000D+00 0.0000D+00 

0. 0000D+00 0.0000D+00 3.0000D+01 

NSCP = 

9. 

DIMS(9) = 

3.1000D-02 
DIMS(IO) = 

5. 

DIMS(11:13) = 

2. 3. 5. 

NMOD = 

1 . 


-3.2163D+01 0.0000D+00 0.0000D+00 

-3.1 000D-1 4 0.0000D+00 0.0000D+00 

-3.9954D-1 3 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 -1.4924D+01 3.3200D+01 

0.0000D+00 -3.3200D+01 -1.4924D+01 
0.0000D+00 0.0000D+00 0.0000D+00 

O.OOOOD+OO 2.1266D+00 -7.9621 D-02 
0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 0.0000D+00 0.0000D+00 

0.0000D+00 O.OOOOD+OO 0.0000D+00 
O.OOOOD+OO O.OOOOD+OO O.OOOOD+OO 
I.OOOOD+OO O.OOOOD+OO O.OOOOD+OO 
2.8682D-1 3 O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO O.OOOOD+OO 


O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
-5.6372D+00 O.OOOOD+OO 
8.2708D+01 O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
I.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO I.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
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LBOUND 


-24. -25. 


UBOUND = 


1.0500D+01 2.5000D+01 
SGDP = 


COLUMNS 1 
9.9858D-01 
-2.2064D-05 
-3.0376D-06 
-4.8550D-08 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
O.OOOOD+OO 
O.OOOOD+OO 
0.0000D+00 
1.0000D+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
6.6066D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


[THRU 6 
-4.9898D-01 
9.9029D-01 
-2.4344 D-02 
-3.7840D-04 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
I.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
2.8002D+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-2.3061 D-02 
3.0419D-02 
9.9263D-01 
3.0888D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
I.OOOOD+OO 
O.OOOOD+OO 
9.7528D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-9.9636D-01 
1.1012D-05 
1.5615D-06 
I.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
I.OOOOD+OO 
2.8682D-1 3 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-6.8684D-04 
-8.951 7D-0 6 
-1.9730D-04 
-2.1 021 D-06 
3.2457D-01 
-5.3951 D-01 
1.1730D-02 
5.8018D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


COLUMNS 7 
-1 .9398 D-02 
-3.5455D-04 
-5.5898D-03 
-1.6246D-04 
O.OOOOD+OO 
O.OOOOD+OO 
-1 .2748 D-01 
-7. 0635 D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
2.651 5D-01 
2.5144D+01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


THRU 12 
-2.3460D-02 
-3.4715D-04 
-6.7448D-03 
-1.1376D-04 
O.OOOOD+OO 
O.OOOOD+OO 
7.0635D-02 
-1 .2748 D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-1.1230D-02 

-2.1587D-04 

-5.1480D-03 

-9.2150D-05 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

3.9455D-01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

5.7684D-02 

O.OOOOD+OO 

3.0000D+01 

O.OOOOD+OO 

O.OOOOD+OO 


1 .4204D-03 
2.2006D-05 
3.0294D-06 
4.8464D-08 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
9.9475D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
I.OOOOD+OO 
O.OOOOD+OO 


1.6756D-03 

3.2595D-05 

8.1752D-05 

1.2718D-06 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

9.9475D-01 

-5.2335D-05 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

3.3670D-03 


-1.9073D-04 
-2.3042D-06 
-5.4769D-05 
-3.9946D-07 
5.3951 D-01 
3.2457D-01 
6.3442D-03 
8.0366D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


4.41 28D-08 
8.4718D-1 0 
2.1547D-09 
2.231 6D-11 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
5.2335D-05 
9.9475D-01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


COLUMNS 13 THRU 16 


-4.6986D-04 

-5.9214D-06 

-1.3495D-04 

-1.2327D-06 

7.7808D-01 

1.8085D+00 

1.2418D-02 

1.5299D-02 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 

0.0000D+00 


-2.0145D-04 
-3.2479 D-06 
-9.2150D-05 
-1.0219D-06 
0.0000D+00 
0.0000D+00 
0.0000D+00 
0.0000D+00 
2.0182D-02 
0.0000D+00 
0.0000D+00 
O.OOOOD+OO 
0.0000D+00 
0.0000D+00 
0.0000D+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


NSGDP = 

12 . 

SDC = 

8.0880D-01 
O.OOOOD+OO 
-8.2059D+00 
-1 .41 07D+00 

NSDC = 

2 . 

YW = 

2. 3. 5. 


-7.8095D-01 

I.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


1.2858D-05 
1 .991 2D-07 
2.8235D-08 
2.9631 D-10 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
1.8013D-02 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


-1.0667D+00 

3.1000D-02 

O.OOOOD+OO 

O.OOOOD+OO 


1.8572D-05 

3.5655D-07 

9.0683D-07 

9.3896D-09 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

2.2037D-02 

-9.3255D-01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


-2.1447D+00 

3.1000D-02 

O.OOOOD+OO 

O.OOOOD+OO 


B.2 TIME EVALUATION 


<> EXEC ('M UGD I R:TIM EAN AL.CM D') 
SELECT COMMANDFILE 

1 TEVAL.CMD 

2 TGUST.CMD 

3 TPLOT.CMD 

4 USERFILE 

5 AJO USERFILE 

6 QUIT 

SELECT> 1 


1.2979D-01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 
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CHOOSE PLANT 

1 CONTINUOUS PLANT 

2 DISCRETE PLANT 

SELECT> 1 
IS LIMITER USED 

1 NO 

2 YES 

SELECT> 2 

SAME SIGNALS FOR ALL MODELS 
SELECT 

1 YES 

2 NO 

SELECT> 1 

PLEASE ENTER TIME MAXIMUM [10] 

100 DATA PT LIMIT OK 

1 YES 

2 NO 

SELECT> 1 

ENTER TIME STEP WHICH IS APPROX. TMAX/100 [DELT] 
SELECT SIGNAL 

1 PULSE 

2 RAMP 

SELECT> 2 

NUMBER OF SIGNAL COMMANDS 
NSIG = 

3. 

SELECT CHANNEL(S) FOR SIGNAL COMMANDS <I1,I2,ETC> [1] 
ENTER RAMP SLOPE(S) <I1,I2,ETC> [.2] 

ENTER RAMP SATURATION <I1,I2,ETC> [1] 


TIME EVALUATION FOR FIRST RUN COMPLETED 


CHANGE CHANNELS 

1 YES 

2 NO 

SELECT> 1 

SAME SIGNALS FOR ALL MODELS 
SELECT 

1 YES 

2 NO 

SELECT> 1 

PLEASE ENTER TIME MAXIMUM [10] 

100 DATA PT LIMIT OK 

1 YES 

2 NO 

SELECT> 1 

ENTER TIME STEP WHICH IS APPROX. TMAX/100 [DELT] 

SELECT SIGNAL 

1 PULSE 

2 RAMP 

SELECT> 1 
DEFAULT MENU: 

1 KEEP DEFAULTS 

2 REVIEW VALUES 

3 START TIME 

4 NON-SELECTED PULSE 

5 PULSE FREQUENCY 

6 SAVE CHANGES 

SELECT> 1 

NUMBER OF SIGNAL COMMANDS 
NSIG = 

3. 

SELECT CHANNEL(S) FOR SIGNAL COMMAND <11 ,12,ETC> [1] 

SELECT PULSE WIDTHS FOR SIGNALS SELECTED <I1,I2,ETC> [10] 

SELECTED PULSE MAGNITUDE(S) <I1,I2,ETC> [100] 

TIME EVALUATION FOR SECOND RUN COMPLETED 
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CHANGE CHANNELS 

1 YES 

2 NO 

SELECT> 2 

DELETE SYSTEM BUILD 

1 YES 

2 NO 

SELECT> 1 

END OF TEVAL.CMD 


SELECT COMMANDFILE 

1 TEVAL.CMD 

2 TGUST.CMD 

3 TPLOT.CMD 

4 USERFILE 

5 AJO USERFILE 

6 QUIT 

SELECT> 3 
PICTURE TAKING: 

1 EVERYTHING 

2 SOMETHINGS 

3 NOTHING 

SELECT> 2 
Available devices are: 

1 VT220 

2 VT240/Monochrome 

3 Tektronix 4014 

4 MATLN03 

SELECT> 4 

ENTER GRAPHICS FILENAME [MATLASER.DAT]: CARRIAGE RETURN 
SELECT ANALYSIS: 

1 TIME EVALUATION 

2 GUSTS 

3 QUIT 


SELECT> 1 


PLOT MENU 

1 WHOLE 

2 1/2 UP/LOW 

3 QUIT 

SELECT> 2 

ENTER SIMULATION NAME(MAX 22 CHAR) FOR PLOT TITLE (NAME SHOULD HAVE AN 
ASTERISK AT THE BEGINNING AND BE SURROUNDED BY APOSTROPHES) 

’‘REGULATOR’ 

NO.COLUMNS.NO. RUNS, NMOD,PLANT OUTPUTS, CONTROLS, MODEL, SIGNAL COMMANDS 
T 

22. 2. 1. 5. 2. 3. 1. 

SELECT FOR UPPER 

ENTER VECTOR OF VARIABLES DESIRED [2 3 411] 

TOP HALF PLOTTED (FIG. 21) 


SELECT FOR LOWER 

ENTER VECTOR OF VARIABLES DESIRED [5 6 7] 

LOWER HALF PLOTTED (FIG. 21) 


TAKE PICTURE NOW: 

1 YES 

2 NO 

SELECT> 1 
SELECT ANALYSIS: 

1 TIME EVALUATION 

2 GUSTS 

3 QUIT 

SELECT> 3 


B.3 TIME RESPONSE TO EXTERNAL DISTURBANCE 


SELECT COMMANDFILE 

1 TEVAL.CMD 

2 TGUST.CMD 

3 TPLOT.CMD 

4 USERFILE 

5 AJO USERFILE 

6 QUIT 

SELECT> 2 
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<Y, T,ST>=GUSTR(SGDP,NSGDP,SDC, NS DC, DELT.YV, COUNTER) 
CHOOSE CHANNELS 


1 CHANNEL1 

2 CHANNEL2 

3 QUIT 

SELECT> 1 
NO. OF DATA PTS. 

1 KEEP DEFAULT (1000) 

2 CHANGE 


SELECT> 1 


TIME RESPONSE FOR FIRST CHANNEL COMPLETED 


CHOOSE CHANNELS 

1 CH ANN ELI 

2 CHANNEL2 

3 QUIT 

SELECT> 2 
NO. OF DATA PTS. 

1 KEEP DEFAULT (1000) 

2 CHANGE 

SELECT> 2 

ENTER NO. OF DATA PTS [2000] 

TIME RESPONSE FOR SECOND CHANNEL COMPLETED 


CHOOSE CHANNELS 

1 CHANNEL1 

2 CHANNEL2 

3 QUIT 

SELECT> 3 

END OF TGUST.CMD 


SELECT COMMANDFILE 

1 TEVAL.CMD 

2 TGUST.CMD 

3 TPLOT.CMD 

4 USERFILE 

5 AJO USERFILE 

6 QUIT 

SELECT> 3 
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PICTURE TAKING: 

1 EVERYTHING 

2 SOMETHINGS 

3 NOTHING 

SELECT> 2 
SELECT ANALYSIS: 

1 TIME EVALUATION 

2 GUSTS 

3 QUIT 

SELECT> 2 
PLOT MENU 

1 WHOLE 

2 1/2 UP/LOW 

3 QUIT 

SELECT> 2 
SIGNAL INPUTS: 

1 CHANNEL 1 

2 CHANNEL 2 

SELECT> 2 
NCY = 

2 . 

NO. COLUMNS.NO. MODELS.TOTAL OUTPUTS, PLANT OUTPUTS, CONTROLS.GUST 
V 


9. 1. 9. 5. 2. 2. 

SELECT FOR UPPER 

ENTER VECTOR OF VARIABLES DESIRED [2 9] 

TOP HALF PLOTTED (FIG. 22) 

SELECT FOR LOWER 
SIGNAL INPUTS: 

1 CHANNEL 1 

2 CHANNEL 2 

SELECT> 2 
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NO.COLUMNS.NO.MODELS, TOTAL OUTPUTS, PLANT OUTPUTS, CONTROLS, GUST 
V 


9. 1. 9. 5. 2. 2. 

ENTER VECTOR OF VARIABLES DESIRED [6] 

LOWER HALF PLOTTED (FIG. 22) 


TAKE PICTURE NOW: 

1 YES 

2 NO 

SELECT> 1 
SELECT ANALYSIS: 

1 TIME EVALUATION 

2 GUSTS 

3 QUIT 

SELECT> 3 

SELECT COMMANDFILE 

1 TEVAL.CMD 

2 TGUST.CMD 

3 TPLOT.CMD 

4 USERFILE 

5 AJO USERFILE 

6 QUIT 

SELECT> 6 


B.4 ADDITIONAL OUTPUT DATA 


GUSTAT = 

COLUMNS 1 THRU 6 


5.2292D-01 6.3624D-04 

-2.9174D-01 2.2937D-04 

5.0874D-04 

-2.0502D-04 

3.2339D-03 
-6. 7829 D-04 

3.2607D-03 

-1.3392D-03 

1.0395D-02 

-3.1300D-03 

COLUMNS 7 THRU 9 
1.7837D-03 1.0000D+00 

-5.3791 D-04 5.9155D-02 

0.0000D+00 

0.0000D+00 




GWSTAT = 





COLUMNS 1 THRU 6 
3.3245D+02 3.8928D-01 

2.3273D+01 -2.2649D-02 

2.8748D-01 

7.8752D-04 

1.0632D+00 

3.0784D-01 

2.3758D+00 

1.0629D-01 

2.0152D+01 

1.4603D+00 

COLUMNS 7 THRU 9 
3.4581 D+00 0.0000D+00 

2.5132D-01 O.OOOOD+OO 

1.0000D+00 

4.3018D-02 




































ANALYSIS 


1. CHANGE OMEG^ 

2. S.V. ERROR 

3. S.V. LOOP TRANSFER 

4. BODE FREQ RESPONSE 

5. SSCOV 

6. CLOSED LOOP EIG 

7. QUIT 


- OMEG MENU 

1. CHANGE OMEG 

2. REVIEW OMEG 

3. QUIT 


Figure 3. Frequency-evaluation menu. 



Figure 4. Multiplicative error at plant input. 



Figure 5. Input-to-output map. 
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FREQUENCY - RAD/SEC 



BODE FREQUENCY RESPONSE 

Figure 18. Computer-generated output for Bode frequency response. 
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Figure 19. Computer-generated output for Nyquist plot. 
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Figure 20. Computer-generated output for steady-state covariance plot. (State and output numbers added for 
clarity.) 
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Figure 21. Computer-generated output for regulator time response. (Numbers on plots represent outputs and 
are added for clarity.) 
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